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Abstract 

A systematic approach to long-distance effects in exclusive radiative weak 
B decays is presented, based on a combination of the heavy quark limit 
with perturbative QCD. The dominant long-distance effects, connected with 
weak annihilation and VF-exchange topologies, can be computed in a model- 
independent way using experimental data on B radiative leptonic decays. 
Nonfactorizable corrections vanish in the chiral limit and to leading twist. 
The left-handed photon amplitudes are shown to be enhanced relative to the 
right-handed ones, both in the long- and short-distance parts of the B de- 
cay amplitudes. Recent CLEO data onB-> K*"f decays are consistent with 
Standard Model estimates of the long-distance contributions, and disfavor an 
enhanced gluonic penguin contribution. We discuss the implications of our 
results for the extraction of \Vtd\- 
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I. INTRODUCTION 



Rare radiative decays of B mesons received considerable theoretical attention due to their 
special sensitivity to physics beyond the Standard Model. Alternatively, if the validity of the 
SM is taken for granted, these decays can offer useful information on the magnitude of the 
CKM parameters. So far, most of the theoretical effort has been concentrated on aspects 
of inclusive decays b — > 57, which can be treated with the help of an operator product 
expansion (OPE) combined with the heavy quark expansion The corresponding strong 
interaction effects can be quantified in terms of a few nonperturbative matrix elements. 
Experimental results on the branching ratio for inclusive B — > X s ^ decays |f2],|3| appear to 
be in agreement with the predictions of the Standard Model. 

On the other hand, the treatment of exclusive radiative decays such as B — > K* / y or 
B — > p7 is considerably more difficult. This is due to the fact that bound state effects must 
be taken into account, which are essentially nonperturbative in nature. In addition to the 
matrix element of the penguin operator (short-distance component), a nonlocal contribution 
must be included too, arising from a product of the electromagnetic coupling with the weak 
nonleptonic Hamiltonian (long-distance component). These effects have been estimated 
using various methods, such as perturbative QCD combined with the quark model |5|,|6|,|9HTT 
dispersion methods T2]| , vector meson dominance [13,13] and QCD sum rules |T5|-|rT||. 



Recently the CLEO collaboration measured the branching ratios for several exclusive 
radiative decays B — > Vj, with V = K*°,K* + ||, and gave upper limits on the V = p, u 
modes. A more detailed study of these decays is therefore well motivated. 

In this paper we propose a different approach to the calculation of the dominant long 
distance contributions to B — > V7 decays (with V a member of the SU(3) octet of vec- 
tor mesons), arising from weak annihilation and W -exchange topologies. We start by 
parametrizing all these decays in terms of a few common amplitudes (for each photon helicity 
X = L,R), using SU(3) symmetry. It can be shown that the structure of these universal am- 
plitudes simplifies very much to the leading order of a twist expansion, in powers of A/E^ 
and A/p + = A/me, with p + the light-cone component of the vector meson momentum. 
Since the photon energy in B — > p(fT*)7 decays is E 1 = 2.6 GeV, the higher twist correc- 
tions can be expected to be well under control. We find that, to leading twist, the photons 
emitted in B decay are predominantly left-handed. Second, the nonfactorizable corrections 
turn out to be computable in terms of the light-cone wavefunctions of the B and V mesons. 
Furthermore, in the chiral limit, the nonfactorizable corrections vanish exactly to one-loop 
order. The remaining factorizable contribution can be determined in a model-independent 
way using experimental data on radiative leptonic B decays B — > ^ev. 

These results can be used to assess the feasibility of the determination of the ratio 
from a comparison of the B — > K*^ and B — > p7 decays. The Appendix describes 
certain constraints on radiative decay form-factors following from a Ward identity expressing 
the conservation of the electromagnetic current. 

II. FLAVOR SU(3) PREDICTIONS FOR B -> V>y DECAY AMPLITUDES 

In the Standard Model, the radiative B — > decays (with V a member of the SU(3) 
octet of vector mesons) are dominated by the operators (q = s,d) (see, e.g., [ |19|1 ) 
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n rad = -^f\[ s) E ci^ol'V) + («-<*), (i) 

V2 i= 7,8 

= r^F^^(m 6 P B + m g P L )6 , (2) 

# = Tll9VK^ + ^)^ & . ^ = ^±75) (3) 



We denote here and in the following the combination of CKM matrix elements Aj^'-* = V q bV* q ,. 
In addition to this, there are also contributions mediated by the usual weak nonleptonic 
Hamiltonian 

AC 

Hw = —^ {\^\C x {ub) v „ A (su) V - A + C 2 (sb) v - A (uu)v- A ] (4) 

10 i 

+ A( s )[C 1 (cfo) v _ A (sc) v _ A + C 2 (sb) v _ A (cc)v_ A ] - A t (s) £ c,gS s) + (s - d) 

1=3 J 

with (gi^V-yife^V-A = [qi7ij,PlQ2)[Q31 >j 'PlC[a]- In these diagrams, the photon attaches 
to the internal quark lines with the usual electromagnetic coupling. We will neglect in the 
following the contributions of the penguin operators Q3-10, which are suppressed relative to 
those of the current- current operators by their smaller Wilson coefficients. 

In a quark diagram language, all B — > V^x decay amplitudes are parametrized by nine 
independent SU(3) amplitudes, for each individual photon helicity A = L, R. They appear 
together in fewer combinations, such that a certain predictive power is still retained. They 
include penguin-type amplitudes P U)C> t - corresponding to topologies with a u-, c- and t- 
quark running in the loop, the weak annihilation (WA)-type amplitude A, the W-exchange 
amplitude E, penguin- annihilation amplitudes PA, and amplitudes with one insertion of 
the gluonic penguin operator C 8 which will be denoted (see Fig. 1). In the penguin 
P®. and AfW amplitudes we distinguish between diagrams with the photon coupling to 
the loop quark or the other emerging light quark (i = 1), and diagrams with the photon 
coupling to the spectator quark (i = 2). We adopt the usual phase conventions for the 
vector states (p + ,p°,p~) = (ud, -^(dd - uu), -du), (K* ,K*~) = (sd,-su), (K* + ,K*°) = 

(us,ds), (j)^ = ^g(2ss — uu — dd), and for the heavy mesons (B + , B°, B s ) = (bu,bd,bs), 
(B-,B°,B 8 ) = (bu,-bd,bs). 

The weak annihilation amplitudes A^ contribute only to the B ± radiative decays, where 
they appear always in the same combination A = \A^ + ^A^ — ^A^ (see Fig. 1(a)). 
Writing explicitly the contributions of the penguin amplitudes one has 

A(B- -> p- lx ) = AW(P« + Q U PS + A x ) + Af(P« + Q U P^) (5) 

+\[ d \p tx + M ( x 1) +Q U M^) 
A{B- -> K*~ lx ) = ^\PS + QuPS + Ax) + X { C S) (P^ + QuP% ] ) (6) 

+A| S) (P tx + M« + Q u Mf) . 

The penguin amplitude with an internal t-quark P t x (arising from the penguin opera- 
tor 7 ) is usually called the short- distance contribution, in contrast to the long-distance 
amplitude, arising from nonlocal insertions of the weak Hamiltonian (H) with one photon 
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attachment. In a hadronic language, the penguin amplitudes P^> receive contributions from 
rescattering effects of the form B ± — ► {p ± p } — > p ± 7, and the charmed penguins PW arise 
from rescattering processes P — > (P* 1 * 1 !)* }, {^p } — > p 7 . 

The P° decay amplitudes contain the W^-exchange topology, parametrized by the graph- 
ical amplitudes (Fig. 1(b)). All the decays discussed below contain these amplitudes in 
the same combination E = |P (2) + |P (3) - |P (1) 

^(5° - p° 7A ) = ^A?(P« + Q d Pi? - P A - (Q w - g d )PA uA ) (7) 

+^AW(P« + Q d P% - (Q u - Q d )PA cX ) + ^\f\Ptx + + Q d M { x 2) ) 

A(B° -> iT° 7A ) = A«(P« + g d Pi?) + A^(P C ( A 1) + Q d P%) + A| s) (P tA + M« + g d Mj[ 2 } ) (8) 
A{ B° - 0( § ) 7A ) = - -L A W (P W + g^pP) + Ea + ( g u + Qd _ 2 g s )PA uA ) (9) 

-4|Af (PW + g d P c ( A 2) + (g M + Q d - 2Q s )PA cX ) - ^t d) (Ptx + + ^ M l 2) ) • 

Finally, the corresponding B s decays are written in terms of the amplitudes previously 
introduced as 

A(B S -> K*° lx ) = -Af(Pii ) + Q S PS) - AW(P C ( A 1} + g s P c ( A 2) ) (10) 
-\f\p tx + M ^ +Q s Mf ] ) 

a(b s - p° 7A ) = -j=A^(p A + (g„ - g d )PA uA ) + -^A( s )(g u - g d )PA cA (11) 

A(B S -> 0( g ) 7A ) = -4|AW (2P« + 2g s Pi 2) - P A - (Q u + Q d - 2Q s )PA uX ) (12) 



--J=AW(2P C ( A 1) + 2g s P c ( A 2) - (Q u + Q d - 2Q s )PA cX ) - \j\>i s \Ptx + m{ 1} + Q s M ( x h 



Note that the long-distance penguin-type amplitudes with an internal u and c-quark are 
different in P and P°, B s decays, due to the different electric charge of the spectator quark 
in the two cases. This is different from the conclusion reached in 0]; present experimental 
data from CLEO (see Eqs. (|56|), (E7I) below) appear to confirm the presence of these long- 
distance contributions at the la level. 

Decays to the SU(3) singlet vector meson = —-^{uu + dd + ss) introduce new 
amplitudes S\, which arise from diagrams similar to PA X (see Fig. 1(d)) but with the 
photon attaching to the quark in the loop or to the spectator in the B meson 

A(B° - 0« 7A ) = -^\ P S + Q* P S + E * + Sux) (13) 

-i=A^)(P« + Q d P% + S cX ) - -^A^(P« + M« + Q d M?) 
A(B S -> 0« 7a ) = 7^Al s )(P« + g d Pi 2) + P A + S uA ) (14) 



+4rf)(p« + g.pi? + s cX ) + -4aJ s) (pW + m! 1 ) + g 



.^/g c V- 1 c A 1 ^a- 1 cA ' ^c-v 1 /g "t V- 1 t\ 1 •"-"■A 



(2). 



The decay amplitudes into the physical states 0, u are obtained by combining the octet 
and singlet amplitudes with a mixing angle Qy as 

= 0( 8 ) C os^-0 (1) sin^ y (15) 
a = (j) {8) sin Qy + (1) cos Qy . (16) 

The physical value of the mixing angle is close to tan Qy = corresponding to ideal mixing 
(f) = ss, to = —-^(uu + dd). Combining the amplitudes for decays into SU(3) eigenstates 
one finds the following amplitudes for decay into physical states u, 4> 

A(B° - u lx ) = -i=AW(P« + Q d P%! + E X + l -PA uX + 2 -S uX ) (17) 
-±=/*\PS + Q d P$ + \PA* + \S CX ) - i=Af (P« + M« + Q.Mf ) 
A(P° - 7A ) = ^(-PAa + Sua) + ^A^(-PA cA + S cA ) , (18) 
and 

A(B S - ^ 7 a) = -^(^a + ^PA uA + ^ uA ) + -La^(^PA cA + ^S cX ) (19) 
A(P S - 7 a) = -\ { :\PS + QsPjS - \PA UX + l -S uX ) (20) 

-Ai s) (^ ( A 1} + QsPS! - \pa cX + \s cX ) - A| S) (P, A + m« + Q s Mf) . 

The radiative decay widths are given, in terms of these amplitudes, by 

r(P^ 7 ) = -^ £ L4(P-.V 7A )| 2 . (21) 



A. The short-distance amplitude 

Due to the smallness of the light quark masses ma, m s appearing in the short- distance 
Hamiltonian ([[]), relative to the 6-quark mass, it couples predominantly to left-handed pho- 
tons. The amplitude for emitting right-hand polarized photons is suppressed relative to the 
one for left-handed photons by the quark mass ratio \PtR./PtL,\ — Ofa^s/mb)- Therefore, 
the right-handed amplitude P t R can be neglected to a good approximation. 

The left-handed short-distance amplitude can be expressed in terms of the B — > V form 
factors of the tensor current, defined as 

(V(p\ e)\qa^b\B(p)} = g + {q 2 )e^ vXa e x {p + p') a + g^(q 2 )e^ Xa e x (p - p') a (22) 

+h(q 2 )e^ Xa (p + p') x (p - p') a (e* ■ p) . 

Using the relation 0^75 = \e iiva fi< jafi (corresponding to e 0123 = 1), one finds 
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(V(p',e)\qa^ lb b\B(p)) = 

~W+{q 2 ) [e*(p + p')u - <£ (p + p%] " ^-(<? 2 ) [e*(p - p% - <(P - P%] 
-ih(q 2 )(e* ■ p)[(p + p%(p - p')„ - (p + j/)„(p - p%] • 



(23) 



The relevant quantity is the formfactor g + at the kinematical point q 2 
computed using light-cone QCD sum rules (LCSR) 



0. This has been 
OJ21| and lattice QCD p2], with results 



which are in reasonable agreement with each other. The most recent LCSR calculation 



gave g { f{0) 



lattice result g^ \o) 



0.29 ± 0.04 and g ( +*\o) = 0.38 ± 0.06 |2T|, which compares well with the 
2*2]] . A somewhat larger value \o) = 0.4 was recently 



n OO+0.04 
u - oz -0.02 



extracted [23[ from the D — > K*eu semileptonic form factors using heavy quark symmetry 



relations J24 



Using the latter value for g+(0), one finds for the short- distance penguin amplitudes with 
an internal t quark 



tL 



1.8 x 10" b ) 



C 7 (m 6 )m () m B i? 7 g + (g 2 = 0) 
m b (GeVY 



167T 2 

6^ f g±B 

0.4 



(24) 



4.2 



GeV 



and P t R — 0. We used in this numerical estimate the leading-log value for the Wilson 
coefficient CV(m&) = —0.31. In the following section we will prove a similar suppression of 
the right-handed amplitude relative to the left-handed one for the weak annihilation (1) and 
the VT-exchange (E) amplitudes in B decays, to leading order in an expansion in powers of 
1/E T 



B. Weak annihilation and V^-exchange amplitudes 

We start by assuming factorization; the corrections to this approximation will be dis- 
cussed below. The factorized weak annihilation amplitude contributing to the B~ — > p~7 
decay is written as 

A A = ^ ai {-/ B p fl >--7A|(du^ , (25) 

corresponding to the photon coupling to the B~ or to the p~ constituent quarks respectively. 
To tree level, the Wilson coefficient a\ is given by a\ = C\{m^) + C2(mb)/N c = 1.02 fli~9"| . 
The first term can be computed exactly in the chiral limit with the result 

p Bt x{p~lx\{du)^_ A \Q) = -(p~<y x \id fl (du)'{ / _ A \0) = -em p -f p -(e* ■ e* x ) +0{m u ,m d ). (26) 

We used here the following relations for the divergence of the isovector and isoaxial currents 
in the presence of an electromagnetic field A (with e > 0) 

id^rf^q) = -eA x q[Q, X]<y x q - ?K A*]g (27) 
id^q^^q) = -eA\q[Q, A*]7a759 - q{™, A*}759 • (28) 
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Q = diag (|, — |, — |) and m = diag (m u ,m d ,m s ) are the quarks' electric charge and mass 



matrix, respectively. The result ( p6[ ) is analogous to a similar one obtained in |25[ for the 
long-distance contribution to the decay B — > 7re + e~. 

The remaining contribution can be expressed in terms of the two form factors fy^E^) 
parametrizing the radiative leptonic decay B~ — > ^ev 

' (l(q,e x )\q lfl (l- l5 )b\B(v))= (29) 



7T« 



v, q)fv{E 1 ) + [e^(w • g) - q p (e* x ■ v)\f A {E 1 ) - - Q b )f B m B {v ■ e* x ) Vfl , 

where v denotes the B meson velocity (ps = vtlbv). The last term is present only for 
charged B mesons, and is required for the gauge invariance of the complete decay amplitude, 
including also the lepton part. Although it is not relevant for real photon processes, it does 
contribute to decays involving virtual photons, which can probe the /i = component of the 
electromagnetic current. A detailed derivation of this term is presented in the Appendix. 

Combining the two terms in (^) one finds the following result for the weak annihilation 
helicity amplitudes A\ 



G 



Al,r = -^ ai em p -f p - [f B + E.U'a W =F f v ° >(E,))\ . (30) 

Eventually the form factors fv,A{E 1 ) will be extracted from the doubly differential spec- 
trum d 2 T / dEedE^ in the radiative leptonic decay B — > 'yeu, which will allow a model- 
independent calculation of the weak annihilation amplitudes Al )B . A considerable simpli- 
fication can be achieved if the WA amplitudes are expanded in powers of K/E 1 using the 



(B-)/ p x _,_ AB-). 



methods of perturbative QCD for exclusive processes ||26|| . It has been shown in [27] that 
the leading terms in the expansion (of order Oil/E^)) of fy^E^) and /^(£^ 7 ) are related, 
and can be expressed in terms of the valence light-cone wave function of the B meson as 

*{B±)(k \ _ J-^J/p \ _ fB m B f n n _Qb 



R is a hadronic parameter given by an integral over the B meson light-cone wave function 
R = J °° dk + ^B{k + ) /k + (with the normalization J °° dk + ^B{k + ) = 1). A similar result is 
obtained in the quark model, with the identification R — > l/m u , with m u ~ 350 MeV 
the constituent light quark mass The relation ( |3~TD among fy^E^) implies that a 

measurement of the photon spectrum in B ± — > 7ez/ should be sufficient for their extraction 
(without any knowledge of the B meson light-cone wavefunction) . 

For the moment, in the absence of such data, the parameter R can be estimated using 
model wavef unctions, which results in typical values around R = 2.5 ± 0.5 GeV -1 |27| . We 
will use this central value in our estimates below. 

It is possible to give a model-independent lower bound for the R parameter, in terms of 
the B meson binding energy A = m# — rrib. At tree- level this bound reads R > 3/ (4A) | |27| |. 
This can be used together with (^) and (|31|) to derive a lower bound on the magnitude of 
the Al amplitude (corresponding to B~ — > p _ 7 decays) 



7 



A L > V2G F a iemp -f p - 3QJ r B = (054 x 10-) (^) (IsJM^D.) ( 3M \ GeV . (3: 
b 9 Jp 8A v ; V1.0/V 175 / \ A (MeV) / v 

We used in this estimate f p - =216 MeV, as determined from the leptonic decay width of 
the p° meson. 

To leading order in 1/Ej, the right-handed amplitude Ar receives contributions only 
from the first term in (130[) 



A R = --L ai em p -f p -f B = -(0.07 x lO" 6 ) GeV, (33) 

where we used the same values for the hadronic parameters as above. Comparing with ( |32|) 
one can see that the left /right-handed WA amplitude ratio is enhanced by a factor of more 
than 8. 

The VT-exchange amplitude E is given by an expression analogous to (^), with a different 
phenomenological factorization coefficient a 2 = ^(mj) + Ci(rrib)/N c = 0.17 [|Tj|]. The result 
corresponding to the decay B° — > p°7 is written as 

£ A = G F a 2 {-/ B pB>%A|(™0^ • (34) 

In the chiral limit only the second term contributes, which gives (with f p o = -^f p ±) 

El,r = -^a 2 em p of p+ E^ff\E 7 ) T fP(E 7 )) . (35) 

The form-factors appearing on the RHS are given by an expression analogous to (^Tj) (the 
sign change is due to our phase convention for the B° state) 

#v 7 ,=-/r^)=-^(Q^-|) +0 (|). ( 36) 

To the leading order in A/E^, one finds the lower bound 

E L > V2G F a 2 em p0fp+ gR = (0.05 x 10 ) [—) I^TS - J [TjMeVj) GeV " 

The ly-exchange amplitude E is suppressed relative to the WA amplitude A by color- 
suppression in the ratio a 2 /a\ and by an additional factor \Qd/Q u \ = 1/2. The right- 
handed Eft amplitude receives contributions only from the higher twist terms in (|36D and is 
therefore suppressed relative to the left-handed amplitude Er by a factor of A/E~, ~ 0.15, 
corresponding to a photon energy in B — > p{K*) r y decays E 1 = 2.6 GeV. This gives the 
estimate E R ~ 0.15£ L > (0.007 x 10~ 6 ) GeV. 

The remaining amplitudes P UyC are considerably more difficult to evaluate. A great deal 
of effort has been put into attempts to calculate them, with methods such as simple vector 
meson dominance 0,[T3|,|rj]], dispersion techniques combined with Regge theory (|T2J and light- 
cone QCD sum rules [p^5| -[T7|. In the following section we will revisit some of these estimates, 
putting them into perspective relative to the effects computed in this section (see Table I). 
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Photon helicity 


\Ptx\ 


\Pcx\ 


\PuX\ 


I^aI 


\Ex\ 


A = L 


1.8 


0.16 


0.03 


0.6 


0.05 


\ = R 





0.04 


0.007 


0.07 


0.007 



Table I. Estimates of the short-distance and long-distance amplitudes in B — > 
fry decays (in units of 10~ 6 GeV). The estimates of the WA and ^-exchange 
amplitudes A x and E\ used R = 2.5 GeV -1 . 



C. Long-distance amplitudes with internal c— and u-quark loops 

The long-distance amplitude P c \ induced by the b — > ccq part of the weak Hamiltonian 
(||) is usually assumed to be dominated by the diagram wherein the photon couples to the 
charm quark loop. In a hadronic language, this contribution arises from the weak decay of 
the B meson into Vi) {n) (with ^ (n) a cc state with the quantum numbers of the photon), 
followed by the annihilation of the tp^ 1 ' state into a photon | i~3yi4 1 . For an arbitrary photon 



virtuality q , this long-distance amplitude can be written as a sum over intermediate states 

q 2 - M 2 + iM n T n 

The lowest-lying state contributing to the sum is the J/ip, with a mass of 3.097 GeV. There- 
fore, for a real photon, the J/ip propagator denominator is large, such that its contribution 
can be expected to be strongly suppressed. These heuristic arguments are supported by an 



explicit QCD sum rule calculation |17|], where the charmed penguin amplitude P c \ is found 
to be less than 5% of the dominant short- distance amplitude PtL- This result is confirmed 
by existing experimental data, as discussed below. 

In the following we will substantiate these qualitative arguments with an estimate of 
the sum fl38|), truncating it to a few lowest states. For definiteness we consider the decay 
B~ — > p~7. The weak decay amplitude in (j38"[) can be estimated using factorization, with 
the result 

A{B- -> p-^ n) ) = -|A«a 2 / V (.)f^ w / V te(e;,p B , Pp ,e;) (39) 



p 



2A 2 (m? 



-K + , m p )A l (mf j)(n) )(e* i) ■ e*) + — Zl^l ( e * . q)( e ^-p B ) \ 
The form-factors appearing here are defined in the usual way (with q = p' — p) 

(p-(p',e)\u^b\B-(p)) = 2V{ ~p te(f,,p,p',e*) (40) 

TUB + TTlp 

(p-(p', e)\u^ l5 b\B-(p)) = 2m p A)(g 2 )^r\ + (m B + m p )A l (q 2 ) ^ - ^r%\ (41) 

-M{q 2 ) 6 . 9 \Pn + P'f,- mB 2 mp ( l l A , 
m B + m p \ * q 2 p / 
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and the decay constant of the ip state is given by {0\c^ p c\i})(q,t)) = m^f^e p . Inserting (p9|) 
into the sum (|38| ) and performing the sum over the ij)^ n ' polarization one finds the following 
result for the long-distance charmed penguin contribution to the B~ — > p~7 amplitude 



P 7aJ 



2m 



B 



y f 

m B + m p ^ ' 



ie(q,e 



Xi V 



,(») (42) 
A 2 (0)[(v.q)(e;.e* x )-(v.e* x )(q.e;) 



To put the result into this (explicitly gauge-invariant) form, one must assume the fol- 
lowing relation between the form factors A\ and A 2 at q 2 
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A 1 (0)=A 2 {0) 



m B - m p 
m B + nir 



(43) 



While the approximate numerical values of these form factors V Bp (0) = 0.34, Af p (0) = 
0.26, A 2 P (0) = 0.22 [^TJ do not satisfy exactly the relation (|43f) , it is worth noting that it 
does hold true to leading order in a large energy expansion in powers of 1/E P . In [[Kj, the 



form-factors for the decay of a B meson into a light pseudoscalar (vector) meson have been 
analyzed using an effective theory for the final state quarks [^TJ. To leading twist, all these 
form factors can be expressed in terms of two universal functions £±(E p ) and (\\{E p ) 



V(q 2 ) 
A2{q 2 ) 



m B J 

-^-C(E f ) 
ttl b + m p 



1 + 



m B J 



m. 



(44) 
(45) 

(46) 



The leading twist form factors satisfy certain relations, which are very similar (but not 
completely identical) to relations previously derived in the quark model |32" |. This is not 



completely unexpected, considering that to leading twist only the meson's valence degrees 
of freedom are relevant, which esentially coincides with the quark model picture. 

Using the relations fl4"5j), ( fjE|) one can see that the condition ( fj3| ) is indeed satisfied, up to 
corrections of higher order in m p /E p , proportional to <V The relations (0), fl46|) imply also 
that only left-handed photons are emitted to leading twist. The relevant helicity amplitudes 
are given by 



_ G F 2m B 

cL,r — Wce^=a 2 ■ 

V2 m B + m p 



£ 7 (F(0)±i 2 (0))E/;,», 



(47) 



Using the realistic values for form-factors quoted above, one finds that the left /right-handed 
helicity amplitudes ratio is enhanced by about a factor of 5. Absolute values for the \P C \\ 
helicity amplitudes obtained by keeping only the lowest two states in ( |38| ) are quoted in 
Table I. The corresponding decay constants can be extracted from their leptonic widths and 
are given by f^w = 395 MeV, = 293 MeV |Lj . 

A similar analysis can be performed for the long-distance amplitude P u \ describing quark 
diagrams with w-quark loops, which is given by a sum analogous to (|38|) 
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A { B~W~Vri=e E Mff^ . (48, 

The weak amplitudes A(5 — ► VV) receive, in general, both color-allowed and color- 
suppressed contributions (in contrast to the decays B — > Vip which are purely color- 
suppressed). For example, in B~ — > p~7 decay, the color-allowed component gives a WA- 
type contribution to the radiative decay. To avoid double counting of amplitudes, we will 
keep therefore only the color-suppressed component of the weak amplitude in ( fiSD (symbol- 
ized by the subscript C). Our analysis is different in this respect from other VMD estimates 
of this amplitude . 



Keeping only the contributions of the lowest two states V = p , u, one finds 

= ^ea^^-E., + A) (V(0) ± *(„)) . (49) 

In writing this result we used the nonet symmetry relation (0\u / ~f fJi u\p ) = (Oluj^ulou) , valid 
in the large N c limit or the quark model. The corresponding numerical results are shown in 
Table I. 

The enhancement of the left-handed photon amplitudes relative to the right-handed ones 
in the long-distance contributions we found is somewhat surprising. Calculations of these 
amplitudes using different methods such as QCD sum rules [pjjj -TTTH an d perturbative QCD 
H find a similar result. Such an enhancement does not necessarily hold in the presence of 
new physics. In certain extensions of the Standard Model, such as the left-right symmetric 
model, the right-handed helicity amplitude can become large (proportional to the virtual 



heavy fermion mass) p3fl . Such effects lead to observable consequences such as significant 



mixing-induced CP asymmetries in weak radiative B — > V'j decays |3J . 

To summarize the results of this section, the rare radiative decays are dominated by the 
amplitudes corresponding to left-handed photons. The dominant long-distance contribution 
to the B~ — > p _ 7 decay comes from the W A mechanism; the corresponding amplitude A\ is 
calculable in a model-independent way from experimental data. The left-handed amplitudes 
satisfy an approximate hierarchy of sizes 

E~P U ~ 0.4P C , P c ~ 0.2A , A ~ 0.3P t . (50) 

The relative magnitude of the components with different CKM coefficients can be obtained 
by combining these estimates for the amplitudes with the corresponding CKM factors. 



D. Nonfactorizable corrections 

The leading corrections to the factorization result (|25| ) for the WA amplitudes come, in 
perturbation theory, from the graphs in Fig. 2, with one gluon exchanged between the B 
and p quarks. 

Let us take for definiteness the final vector meson to be moving along the positive z 
axis. Momentum conservation ttibv = q + p is expressed in terms of light-cone components 
as p + = in b and mg = 2E~ f + p_ (we define the light-cone coordinates as p± = p° ± p 3 ) . 
Then to leading order 0(1) in an expansion in powers ofA/p + = A/ mg, the p meson can be 
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represented by its valence component \q{pi)q{j>2)) with parallel momenta p\ = (yp+, 0_, 0j_) 
and p 2 = ((1 — y)p+, 0_, 0j_). In a tensor language this corresponds to the p wave function 



p)a/3 



8V% 



p+(t- f) a p<i>±{y) 



(51) 



where 4>±{y) is the twist-2 chiral-odd structure function as appropriate for a transversely 
polarized p meson. It satisfies the normalization condition 



dy(f>±(y) 




where the decay constant fj = (160 ± 10) MeV p8| is defined as 



(Olua^dlp (p, e)) = ifj(e^p u - e v p^) 



The B meson is represented by the tensor wavefunction [27 



V2k 



' Xc ip B (k + , k ± ) I (k+ + a ± ■ k±)A + ^- 



-75 



A. 



a/9 



7-7+ 

4 



(52) 



(53) 



(54) 



In contrast to the light mesons tt or p, the leading twist description of a heavy meson involves 
in general (beyond tree level) also the transverse motion of the heavy quark. 

Each of the four diagrams in Fig. 2a-d contains IR divergences; however, an explicit 
calculation shows that they cancel in the sum. Therefore, (at least to one-loop order) the 
leading-twist nonfactorizable contribution is computable in perturbation theory as a convo- 
lution of the B and p light-cone wave functions with a IR-finite hard scattering amplitude 
T H 



roo rl 

Sn.fAL = / dk + / dy4>±(y, p)ip B (k + , p)T H (k + , y, p) + higher twist . 
Jo Jo 



(55) 



(This is analogous to a result recently found in []35| for nonleptonic decays B — > Tin.) 
However, it is easy to see that the hard-scattering amplitude vanishes in the limit of massless 
final quarks Th = 0(m u ^). Technically, in this limit, the p side of the matrix element 
contains only traces of an odd number of 7 matrices, which vanish. 

None of these conclusions depends on the fact that the photon is on-shell q 2 = 0, such 
that a similar result can be derived for B — > p±e + e~ decays, mediated by a virtual photon. 
However, a nontrivial nonfactorizable correction is obtained for a longitudinally polarized p 
or a pion in the final state. 

In conclusion, the one-loop nonfactorizable correction vanishes identically in the chiral 
limit and to leading-twist order. Therefore the factorized result (|2"5] ) can be expected to give 
a very good approximation to the weak annihilation amplitude A. 



III. APPLICATIONS 
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A. Implications from experimental data 



The CLEO Collaboration [Q recently reported the branching ratios for the B — > K*^ 
exclusive modes 

B(B ± K*^) = (3.76i°;| ± 0.28) x 10" 5 (56) 
B(B° -> A^ 7 ) = (4.55ig;£j ± 0.34) x 10~ 5 . (57) 

These results give the following ratio for the charge-averaged radiative decay widths 



1.29 ±0.55. (58) 



We used here the lifetime ratio of charged and neutral B mesons t(B ± )/t(B°) = 1.066 ± 
0.024 |6 . 



Isospin symmetry constrains the short-distance amplitude P t i to be the same in both 



decays (|56"D , fl5"7"|). Therefore the deviation of the ratio flSSf) from 1 can only come from the 
long-distance contribution M^ 2 ) — Pi 2 ) arising from the photon coupling to the spectator 
quark. Using the unitarity of the CKM matrix one can rewrite = — — \[ in ([]) and 
(§). Furthermore, neglecting the small CKM factor A|f\ these amplitudes can be written to 
a good approximation as 

A{B~ -> K*- 7 ) = A^ } (p t + (M« - Pi 1 )) + -(M< 2 ) - Pi 2 ))) (59) 



A(B° -> ir*° 7 ) = A| s) (p + (M« - Pi 1 )) - i(M( 2 ) - Pi 2 ))) . (60) 

The short-distance amplitude Pn alone would give the following branching ratios 

B s . d .(B ± -> iT^) = 4.6 x 10" 5 , B sA .(B° -> ir*° 7 ) = 4.5 x 10~ 5 , (61) 
where we used the estimate for P t L shown in Table I together with A^ = 0.039. The B 



meson lifetimes were taken as t b ± = 1.64 ±0.02 ps, r B o = 1.55 ±0.03 ps |36fl. Assuming that 



the central values of the experimental numbers (0), Q5"T ) will be confirmed by future, more 
precise determinations, one is led to conclude therefore that the long-distance effects interfere 
among themselves destructively in B° decays, whereas in B ± they could be significant, and 
contribute up to 10% of the total amplitude. This agrees with the theoretical estimates of 
and disfavor a non-SM explanation of the isospin breaking in the ratio (BBf) based on an 



enhanced gluonic penguin as proposed in [37 



In the SU(3) limit the long-distance contributions to B° and B s decays are the same 
(see flT|)-(|l2T)). Therefore one can use the observed B° — > K* 0, y branching ratio to make 
predictions for other strangeness-changing radiative decays of these mesons. Neglecting the 
OZI-suppressed amplitudes S c \ and PA C \, one obtains the following prediction (together 
with r Bs = 1.46 ± 0.06 ps pi) 



B(B S -> 7 ) = —B(B° -> K*° 7 ) ~ 4.8 x 10~ 5 . (62) 
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Using the estimate for El in Table I, and neglecting again the contributions of the OZI- 
suppressed annihilation penguin PA and SU(3) singlet S amplitudes, one can predict very- 
small rates for the B s modes (we used here \V u t,\ = 1.6 x 10 -2 ) 

B(B S -> p° 7 ) ~ B{B S -> cj 7 ) ~ 0.3 x 10" 8 . (63) 



B. Determining \Vtd\ 



p ± 7 is 



The main interest in observing the exclusive charmless radiative decay £P 
connected with the possibility of determining the CKM matrix element \Vtd\- The amplitude 
for this decay is related by [/-spin symmetry to that for the b —>■ s'j mode B ± — ► K* ±r y 



A(B~ - p~ lL ) = A<V + XfpL 
In the SU(3) symmetry limit cll = cl'l 



A(B- -> K*~ lL ) = X[ s 



*L + K 'v L 



(64) 



A L - P c l, Pl = p' L = PtL + M L - P cL . In the absence 

WaWi = (up to 



of the a, a' terms, the ratio of the two amplitudes is equal to \\ t /X 
SU(3) breaking corrections), which offers a way for determining \Vtd\. 
In the general case, however, the (p~7l) amplitude can be written as 



A(B~ p- lL ) 



\K 



(65) 



with ea^ a = a/p — Ai/PtL (where, following the estimates of Sec. II we neglected the 
charmed penguin P c l and gluonic penguin Ml amplitudes relative to the short- distance 
amplitude P t L and the WA amplitude Al). The estimates of the preceding section give 
Ea — 0.3, and <$a = to leading twist. Furthermore, global analyses of the unitarity triangle 



suggest that the ratio of CKM factors appearing in this amplitude is subunitary [38 



|Al d) | 



V, 



ub 



V, 



cb 



V td 



0.4 



(66) 



The estimates of Sec. II. B show that the total rate for B~ (B + ) radiative decay is 
dominated by left-hand (right-hand) polarized photons. Neglecting the small right-handed 
(left-handed) component, which gives only a second order contribution to the ratio, one finds 
for the ratio of the charge-averaged rates 



B(B- -> p-j) + B(B + -> p+ 7 ) 
B(B- -»• AT*- 7 ) + B(B+ -> AT*+ 7 ) 



(67) 



\p\ \Vt 




\\{d)\ (\\{d)\ N 

! W\ £a cos a cos ^ A + \ 7^dh £A 



\p'\ \Vts\j | \x { ; 

The factor in the braces is bounded from above and below by using the inequality 



l+x + x 2 <l — x cos a cos 4>a + x 2 < 1 — x + x 2 



x 



\X) 



(d) 



-e a ^ 0.12. 



(68) 
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Assuming, conservatively, x = 0.2 (which is on the upper side of the estimates for this 
parameter) gives an uncertainty of about 20% in the determination of \V t d\ arising from the 
factor in braces in (|67|) . The SU(3) breaking factor in the ratio of the penguin amplitudes 
Ipl/lp'l — g+ (0)/<7+ ^(0) = 0.76 ± 0.22 is known with a uncertainty of about 30% (we used 
here the results of the LCSR calculation [plfl). Future lattice QCD calculations will hopefully 
improve the precision with which this ratio is known. Adding these errors in quadrature one 
finds a theoretical error in the determination of \Vtd\ from the ratio (|67D of about 35%. The 
leading twist result <$a = implies that the CP asymmetry in B ± — > p ±r y can be expected 
to be very small (A C p oc sin^)- 

The above estimate for the long-distance contribution ea is in general agreement with 
the quark model of |9|,|T4l and QCD sum rules |T5| , p!6| calculations. The advantage of our 
approach will become apparent once experimental data on radiative leptonic B decays be- 
come available. Such data could be used, as discussed in Sec. II, to actually determine the 
WA amplitude A L , and thereby the long-distance contamination ea- 



IV. CONCLUSIONS 

Present experimental data from CLEO on exclusive radiative weak decays |4| are suffi- 
ciently precise such that the long-distance effects are beginning to be observable. A good 
control over the magnitude of these effects is important for an assessment of the uncer- 
tainty they induce into the extraction of \Vtd\ by combining B~ — > p~7 and B~ — > K*~^y 



decays P,|14fl. Unfortunately, these effects have proved notoriously difficult to treat in any 
systematic way. 

In the present paper we focus on a certain class of long-distance contributions, arising 
from weak annihilation and W-exchange quark diagram topologies. We argue that the 
former gives the dominant long-distance correction to the amplitude for B~ — > p~7 decays. 
These corrections can be computed reliably using factorization, in terms of form-factors 
observable in radiative leptonic decays B + — > / ye + u. The nonfactorizable corrections are 
shown to be very small, as they appear only at nonleading twist. A similar method can 
be used to compute the VF-exchange-type long-distance contribution, relevant for the weak 
radiative B d and B s decays. 

Furthermore, to the leading order of an expansion in powers of A/£" 7 , one can show 
that the coupling of left-handed photons in the long-distance amplitude is greatly enhanced 
relative to that of right-handed photons, just as in the short-distance part of the ampli- 
tude. Such an enhancement had been previously observed in QCD sum rule calculations of 
the long-distance amplitude |15| |; our approach clarifies the theoretical limit in which this 



enhancement holds and quantifies the magnitude of the corrections to it. 

Our results should allow one to reduce the model-dependence of the leading long-distance 
effects in B~ — > p~7 decays, and achieve a better control over the theoretical error in the 
corresponding determination of \V td \. 
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APPENDIX A: WARD IDENTITIES FOR LONG-DISTANCE MATRIX 

ELEMENTS 



We present in this Appendix a set of constraints on certain long-distance contributions 
to weak radiative decays of B mesons. These constraints follow from a Ward identity 
expressing the conservation of the electromagnetic current. Let us consider the following 
matrix element of a local operator O, to first order in electromagnetism and to all orders in 
the strong coupling 

Mq,e)f\O{0)\B{v)) = -ie% J d 4 xe^(f\Tf^{x) O(0)\B(v)) , (Al) 

for any hadronic final state /. The operator O can be any quark bilinear including gluon 
fields, or even a four-quark operator. The electromagnetic current includes contributions 
from both the light and heavy quarks j^ m ' = qOpf^q + f^T^c — ^bj^b. A relation similar to 
( |Al"| ) can be written for the matrix element with the photon replaced with a dilepton pair 
in the final state, coupling through a virtual photon to the hadronic system. 

The conservation of the electromagnetic current implies, in the standard way, a Ward 
identity for the matrix element of the time-ordered product in ( |A1| ) 

-iq,j d 4 xe i ^{f\Tj^(x)O(0)\B(v)) = J d 3 xe~^(f\ , 0(0)] \B{v)) . (A2) 

The commutator on the RHS is nonvanishing only if the operator O carries an electric 
charge, as in the case of decays induced by the weak charged current such as B + — ► 7e + z/, 
B — > 7r(p)7e + z/ or B — > D (y *^e + u. In the following we analyze the simplest such case, the 
radiative leptonic decays of a B meson, for which O = 67^75? and |/) = |0). 

The equal-time commutator on the RHS of (|A2|) can be computed explicitly, with the 
result 



- tq,J d*xe***(0\Tj?*(x) (b lul5 q)(0)\B(v)) = {Q b - Q q )(0\b^y 5 q\B{v)) (A3) 

= {Qb - Qq)fBm B V^ . 

The most general parametrization of the matrix element on the LHS can be written in terms 
of five form-factors fi{q 2 , v ■ q) 



d 4 xe^(0|Tj- m -(x) {ln^ 5 q){0)\B{v)) = + f 2 v„v v + hq,q u + hq,v u + fcv,q v . (A4) 

The Ward identity ( |A^| ) implies two constraints on these form-factors 

{v ■ q)h + q 2 U = {Qb ~ Q q )f B m B , h + q 2 f 3 + {v ■ q)f 5 = . (A5) 

For the case of a real photon q 2 = 0, these constraints fix uniquely the form-factor f 2 {0,v-q), 
and relate /i(0, v ■ q) and /s(0, v • q). From Eq. ([Al] ) one finds thus the following result for 
the matrix element of the axial weak current 

(7(9, e)f\h^q\B{v)) = -f B [(v ■ q)e; - {v ■ e*)q»] + (v ■ e*)v^{Q b - Q q )f B m B , (A6) 
which is the same as the result presented in text Eq. (|29|) , with the identification /s = Ja- 
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FIGURES 

b q 3 b q 3 




qi q2 qi q2 



(b) 




* X * 

(e) 

FIG. 1. Quark diagrams contributing to B — > decays. The cross marks the attachment of 
the photon line, a) weak annihilation amplitudes with i = 1,2,3 corresponding to the photon 
attaching to each of the three quark lines qi,q2,Q3', b) IF-exchange amplitudes E^ l \ i = 1,2,3; 
c) penguin amplitudes Pjp (the photon is attached to the q = d,s quark or the quark q' = u,c 

(2) 

running in the loop) and P\ with the photon attached to the spectator quark; d) annihilation 
penguin amplitudes PA„t\ e) amplitudes with one insertion of the gluonic penguin MM (photon 
attaching to the q line) and (photon attaches to the spectator quark). The box denotes one 
insertion of the operator Qg. 
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(c) 




FIG. 2. Corrections to factorization in B — > V7 decays. Only those diagrams are shown 
which contain IR divergences. The total IR divergence cancels in their sum. The cross denotes the 
attachment of the photon line. 
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